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Abstract. In this paper, we give a unified construction of the recursion operators from the Lax 
representation for three integrable hierarchies: Kadomtsev-Petviashvili (KP), modified Kadomtsev- 
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Petviashvili (mKP) and Harry-Dym under n-reduction. This shows a new inherent relationship between 
them. To illustrate our construction, the recursion operator are calculated explicitly for 2-reduction 
O ' and 3-reduction. 
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1. Introduction 

The recursion operator <£, firstly presented by P.J. Olver pQ, plays a key role (see [2H1] and references 
■ therein) in the study of the integrable system. For single integrable evolution equation, it always 
owns infinitely many commuting symmetries and bi-Hamiltonian structures [2H1] which the recursion 
operator can link. As for an integrable hierarchy, the higher flows can be generated from the lower flow 
. with the help of the recursion operator, which offers a natural way to construct the whole integrable 
hierarchy from a single seed system (see [2H1] and references therein). By now, much work has been 
done on the recursion operator. For example, the construction of the recursion for a given integrable 
system [5rll7|. and the properties of the recursion operator [18H23J. In general, the recursion operator 
has non-local term. So it is a highly non-trivial problem to understand the locality of higher order 
symmetries and higher order flows generated by recursion operator [2?,, 24j . In this paper, we shall 
focus on the construction of the recursion operator and explain the locality of their higher flows 
although the recursion operator associated is non-local. 

The main object that we will investigate is three interesting integrable hierarchies, i.e. Kadomtsev- 
Petviashvili (KP), modified Kadomtsev-Petviashvili (mKP) and Harry-Dym hierarchies [25,26j, which 
are corresponding to the decompositions of the algebra g of pseudo-differential operators 

g:={J2 *H&] = {J2 u i di } © {J>^l := 9>k g <k (1) 

i<SCoo i>k i<k 



X 
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for k = 0, 1, 2 respectively, where Ui are the functions of t = (t\ = x,tz,- • •) and d = d x . The algebraic 
multiplication of d" 1 with the multiplication operator u are defined by 

d i u = ^C J i u {j) d i - j , i£Z, (2) 
i>o 

where u^-* = with 

i(i — 1) ■ - ■ (i — i + 1) 



In fact, 5> fc and are the sub-Lie algebra of 5: #>fc] C #>fc and [g<k,9<k] C g<fc when 

A; = 0, 1, 2. The projections of L = ^ Ui<9* G g to g>fc and are 

L> k = ^Uid\ L <k = ^Uid\ (3) 

i>fc i<k 



Then according to the famous Adler-Kostant-Symes scheme [27J, the following commuting Lax equa- 
tions [25j|26] on g can be constructed 

L tm = [(L m )>k,L], (4) 

where k = 0,1,2 are corresponding to KP, mKP and Harry-Dym hierarchies respectively, with the 
Lax operator L given by 

d + U2&- 1 + u 3 d~ 2 + ■■■ k = 0, 
L=< d + ui + u 2 d~ l + ■■■ k = l, (5) 
uod + u\ + U2&~ 1 + ■ ■ ■ k = 2. 
For simplicity, we rewrite ([5]) in a unified form |25p26| 



i.e. 



l = J^a 1 -', (6) 



l>0 



1-k 



L = Y,u l d 1 - l + nd l ~\ (7) 



l>2-k 



and let 



B m = (L m )>k, L m =Y, = E (8) 

j<m j<m 

Then (j4|) becomes into 

L tm = [S m ,L]. (9) 

These three integrable hierarchies have been studied intensively in literatures [26,28-30 , which contain 
the following well-known 2+1 dimensional equations 

k = : 4u 2 tx = (U2xxx + 12U2U2x)x + 3u 2yy , (KP) 

k = 1 : 4u ltx = (ui xxx - 6uju lx ) x + 3u lyy + 6ui x ui y + Guixxd^uiy, (mKP) 

k = 2: 4u 0t = ulu 0xxx - 3—(uld~ 1 ( — )y) y , (Harry-Dym) 

no «o 
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where we have set £2 = V, *3 = t- 

There are some inherent relationships discovered among these three integrable hierarchies. For 
example, their flow equations are defined by a unified Lax equation ([9]) although their B m are different, 
their Hamiltonian structure is given by a general way, i.e. r-matrix method [28], and there exists an 
interesting link among them in the flow equations and gauge transformations [26j . So it is very natural 
to ask whether there exists a unified way to deal with their recursion operators, which is just our 
central aim of this paper. For the KP hierarchy, W.Strampp & W.Oevel [9] and V.V.Sokolov et al [11] 
separately developed a general method to construct the recursion operator by the Lax representation 
©. V.V.Sokolov et al |11H14] used an important ansatz B = VB n + R that relates B n operator for 
different n, where V is some operator that commutes with the L operator and R is the remainder. 
While, W.Strampp & W.Oevel derived a general expression (see eq.(47) of reference [U]) for the 
recursion operators of the KP hierarchy under n-reduction starting from Lax equations. In this paper, 
we will use W.Strampp & W.Oevel's method. However, their method is not applicable to get a similar 
and compact formula for the mKP and Harry-Dym hierarchies due to following two observations: 



1) (L m ) <0 = J2j<o djb j( m ) for the KP hierarchy, but (L m )<i = J2j<i djb j( m ) + Ejli b j \ m ) for 
the mKP hierarchy, (L m )<2 = £j < 2^"'^i( 771 ) + £^=2 ( m ) + (jn)d\ for the Harry-Dym 



hierarchy. 2) It is not affirmative to get a compact form of the flow equations of mKP hierarchy and 
Harry-Dym hierarchy as eq.(6) and eq.(17) of reference [9] for the KP hierarchy because of the second 
summation terms in the last two cases of 1). 

In this paper, to further find inherent relations between above three hierarchies, we shall improve 
W.Strampp & W.Oevel method (use aj(m) only ) and give a unified construction of the recursion 
operators from the Lax representation for three integrable hierarchies: KP, mKP and Harry-Dym 
under n-reduction (see eq. (fl"8|) ). There are two advantages in our construction: 1) it is easy to 
explain why nonlocal recursion operators produce local flows, since the L.H.S. of © only produces 
the differential polynomials of Ui, thus the flow equations of ([9]) are naturally local; 2) a formula of 
the recursion operator for arbitrary n-reduction are derived, which shows the existence of recursion 
operators for the three kinds of integrable hierarchies, and provides a constructive way to get recursion 
operators for higher order reductions although the calculation is not an easy task. 

This paper is organized as follows. In Section 2, we rewrite the unified Lax equations ([9]) into matrix 
forms in terms of aj(m) under n-reduction. Then, we devote Section 3 to deriving the formulas of the 
recursion operators for the three integrable hierarchies. At last, we consider the applications of the 
formulas of the recursion operators and check the correctness of the formulas. 
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2. Lax Equations 

In this section, we will rewrite the Lax equations ([9]) into matrix forms in terms of dj(m) under 
n-reduction. For this, we start from the mth power of L, that is 

L m = Y,aj(m)d j . (10) 

Thus 

771 

(L m )>k = X>( m )^'' ( n ) 

j=k 

(L m )<k = ^(m)^'. ( 12 ) 

j<k 

Note that, the Lax dynamics equation ([4]) can be rewritten into 

L tm = [L,(L m )<k]. (13) 
We first derive the flow equations for the coordinates Uj. After inserting (|12p into (|13p . we find 

L(L m ) <k -(L m ) <k L = ^^(^-'ajH^'-^-M^a 1 -') 

2>o i<fc 

= E E E(^°fM " CjayM^)^^ 

2>0 i<fc p>0 

= E E E( c rM 9 ~V) - cf^Mu^^ 1 -^- 

g>0 i=0 jr'<fc 

k—1 j+r 

= E EEiW^h-^h^" 

r>l— fc j=—r 1=0 
r>2-k j=l-k 1=0 

According to 0, we know 

^ = E "i,*-^"'- ( 14 ) 

«>2-fc 

So by comparing (fH|) with L(L m ) < / c — {L m ) <k L, we obtain 

r 

u r,t m = ^ O r ja-j(m), r = 2 - fc,3 - fc, • • • , (15) 
j=l-k 

with O rj - given by 

Orj = Y^KT^-i- 1 - C r _- j - l 4- j - l) ). (16) 
1=0 

In particular, we find 

O nr = 0, O r , r -i = u 3 - (1 - r)u 0x . 
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(17) 



Notice that from f|10|) . we find aj(m) can be uniquely determined by Ui, that is, 

a s (m) = mu™~ l u m - s + f S m(u , « m _ s _i), 

where / sm are the differential polynomials in Mo,...,M m _ s _i. After substituting (fTT|) into (fl"5j) . we 
obtains a series of evolution equations for itj. These flow equations are all local because a s (m) are the 
differential polynomials of it,. 

We next consider the so-called n-reduction, that is, we impose the constraints below on the Lax 
operator L , 

L n = (L n )> k . (18) 

Under the constraints above, a s (n) = for s < k. Hence from ()17p . we can express Uj for j > n — k in 
terms of (u2-fc, • • • , u n -k)- Thus only n — 1 coordinates (it2_fc, • • • , u n -k) are independent, 

which are in one-to-one correspond with (afc(n), afc + i(n), • • • ,afc +n _2(n)). For example, under the 2- 
reduction, only «2-fc is independent, then the flow equation (|15p implies the following 1 + 1 dimensional 
equations, 



for k = 



u 2t 5 



for = 1 



"Its 



1 

-^zxx + 3U2-U2X, 

15 2 5 5 1 

— U 2 U2x + ~tU2U2xxx + -^2xx^2x + — ^rxxxx, 

2 4 2 lb 



1 3 2 

15 „ 5 1 



-u x u lx 



^lxxx^l ~i~ n Uixxxxx 

lb 



L lx ~ ~^ u l u lxxU\xi 



for 



U()t, 



u ot 5 



^UqUq xxx , 



(19) 
(20) 

(21) 
(22) 

(23) 
(24) 



After the preparation above, under n-reduction we can at last rewrite the Lax equations ([9]) into 
matrix forms in terms of aj{m). For this, we denote 

U(n) = (u 2 -k,u 3 - k ,- ■ ■ ,Un-kY, 
A(n,m) = (a_ 1+fe (m),o_2 + fc(m),-- 

/ 02-Jfc,l-fc 



0(n) 



\ 




(25) 
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On— k,n— 1— k J 



where t denotes the transpose of the matrix, then we can rewrite (|15p into 



U{n) tm = 0(n)A(n,m). (26) 
It is trivial to know that all of the flow equations in U(n)t m are local, including those in U{n)t m+jn - 

3. Recursion Formulas 

In this section, we will construct the recursion operator. To do this, we have to first obtain a 
recursion formula relating A(n, m) and A(n, m + n) under n-reduction constraint, that is, we try to 
seek an operator R(n), s.t. A(n, m + n) = R(n)A(n, m). 

For this, we consider the relation L m+n = L m L n = L n L m . Assuming n-reduction, we find 



{L m L r 



Comparing with 



we find 



with 



)<k = {^2 ^ a A m )^ a li n )^)<k 
j<k-l l=k 
n 

= ( E EE^H«i W (n^<» 

j<k-l l=k p>0 

n 

= ( E E E C^ a] (m)a?-«\n)d^)<k 

j<k—l q<n l=max(k,q) 
n n 

= E E E cj^woj'-"^. 

j<k— 1 q=j+l— k l=max(k,q) 



j<k-i 

n 

aj(m + n) = Pj q (n)aj- q (m), j < k - 1, 

q=j+l-k 



pj q n= E c i -yr q \n), j<k-i, 

l=max(k,q) 

In particular, we have 



Pjn(n) = Uq, Pj,n-i = On-i(n) + ( j - n + l)a n {n) x . 
We next introduce the (n — 1) x (n — l)-matrix S(n) and the (n — 1) x n-matrix T(n) 
/ Pk-i,o(n) Pk-l,l( n ) F fc-l,n-2(^) \ 



S(n) 



Pfe-2,-i(n) P fc _2,o(n) 



\ -P-n+l+fc,-n+2(w) P- n +i+fc ) _ n +3(n) 
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Pfe-2,n-3(^) 



P-n+l+fc,o( n ) / 



(27) 



(28) 



(29) 



j + l- fc,j + 2- ,n. (30) 



(31) 



(32) 



T(n) 



( Pfe-i,n-i(n) P fc _i >n (n) 

Pk-2,n-2{n) Pk-2,n~l{n) Pk-2,n( n ) 








(33) 



\ P-n+l+fc,l( n ) P-n+l+fc,2( n ) P-n+l+fc,3( n ) "' P-ri+l+fc,n ( n ) / 

So ([221) for j = — 1 + fc, — 2 + • • • , — n + 1 + k can be written into 

A(n, m + n) = S(n)A(n, m) + T(n)(a_ n+fc (m), a_ n _i +fc (m), • • • , a_ 2 n+i+fc(m))*. (34) 

So now the only thing that we need to do is to express (a- n+ k(m),a^ n ^i + k(m), • • • , a_2n+i+fc(^))* 
in terms of A(n,m). For this, we will use the relation L m L n = L n L m . 

n 

(L n L m )<k = £E atW&ajW&Uk 
l=k j<k-l 

n I 

= (EE 52c; ai {n)aM(m)8i«-'). 

l=k j<k-l s=0 

n n—n 

= (EE E C s s+ ^a s+ ,{n)af{m)di + n<k 

j<k—l fi=0 s=max(k~ii,0) 
n n—fl 

= E E E w^m^ 

j<k— 1 fi=0 s=rnax(k—[M,0) 



)<k 



Comparing with (|28p . we obtain 



j(m + n) = E^( n ) fl j-/i( m )i J<fc-1, 



with 



n—/j, 



Q M (n)= E ^ +/ ,a s+ »d s , < // < 

s=mai(f:-/i,0) 

In particular, 

Q n (w) = Uq, Q„_i(n) = na n (n)d + o n _i(n). 
Thus using (f29|) = ([35]) . we obtain 

/ a_i +fe (m + ra) \ 
a-2+k(m + n) 



n. 



(35) 



(36) 



(37) 



a-n+i+fc(™ + n) 
y a- n+k (m + n) J 
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( P-i+k,o(n) 
P-2+k,-i(n) 



-P-i+fc,i( n ) 

P-2+k,o( n ) 



■i+k,n-3(n) P-i+k, n -2(n) \ ( a-i+kim) ^ 
„., P „., „(<r,\ a-2+kijrn) 



+ 



-P-n+l+fc,-n+2( n ) -P-n+l+fc,-n+3( n ) 
P-n+k,-n+l{n) P-n+k,-n+2(n) 

P-l+k,n-l(n) 
P-2+k,n-2(n) P-2+k,n-l(n) ■■■ 



P-2+k,n-4:(n) P~2+k,n-3( n ) 

P-n+i+k,-i(n) P~n+i+kfl{n) 
P- n +k,-2{n) P- n+k -i{n) J 

\ / 





P-n+l+k,l(n) P- n+ i +k ,2( n ) "' P-n+l+k,n-l(n) 







\ P -n+k,o( n ) P-n+k,l{n 



/ 



+ 



/ a_ n _i +fc (m) \ 
a-n-2+k(m) 

a-2n+i+k(m) 
\ a_ 2n+fe (m) y 

na n (n)d + a n _i(n) 

Qn-2{n) 



Qi(n) 
Qo(n) 

a- n+k {m) \ 
a_ n _i +fe (m) 



/ Qo(n) Qi(n) 
Q (n) 



P-n+k,n-2(n) P-n+k,n-l( n ) j 
Qn-3{n) Q n -2(n) \ ( 

Q n -4(n) Q n _ 3 (n) 



a-n+2+k(m) 
\ a_ n+ i +fc (m) J 

a_ n+k (m) \ 
a- n - 1+k (m) 


















Qo{n) 




a~2n+2+k{m) 

\ a-2n+l+fc("l) / 

a_i +fe (m) \ 
a_ 2 +fe(m) 







a-2n+2+A:("l) 
a-2n+l+fe("l) / 



na n (n)d + a n _i(n) 



Q 2 (n) 
Qi(n) 

/ a_ n _i +fc (m) \ 
a_„_ 2+ fc("i) 

a_ 2 „+l+fc(m) 

^ a-2n+k{m) y 








a_ n+2 +fc( m ) 

y a_ n+ i +fc (m) y 






na n (n)d + a n _i(n) 

Qn-2{n) na n (n)d + a n _i(n) y 



So if we denote M(n) and iV(n) asnxn and nx(n-l) respectively, that is 



M (n) = 

/ -na n (n)d - (n - k)a n (n) a 

D- 2 +k,n-2{n) 







-na n (n)d - (1 - k + n)a n (n) a 








\ D- n+kfi (n) D- n+k ,i(n) ■■■ -na n (n)d - (2n - 1 - k)a n {n) x J 
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/ D_i +fc>0 (n) 
P-2+k-i{n) 



N(n) 



D-l+k,i(n) 
D-2+k,o(n) 



P-n+l+k-l( n ) D-n+l+kfi( n ) 

P_ n+ fe_ 2 (n) P_ n+Aj _i(n) y 



P-n+l+fc,-n+2( n ) P-rt+l+fe,-n+3( n ) 
\ P-n+k,-n+l(n) P_ n+ fe_ n+2 (n) 

with Dj S = Pj )S {n) — Q s (n), then we have 

M(n)(a_ n+fc (m), a_ n _i +fc (m), • • • , a_ 2n+2+fc (m), a~2n+i+fcO))* 
= -iV(n)(a_i +fc (m),a_2+fc(m), • • • , a_ n+2+fc (m), a_ n +i +fc (m))*. 



(38) 



Since M(n) is invertible, we can solve (a_ n+ fc(m), a_ n _i + fc(m), • • • , o_ 2n+2+ fc(m), a_ 2n _|_i + &(m)) from 
(1381) and then insert into (OH 



So we get 



with 



If we set 



then we can easily find 



U(n) 



■m-j-jn 



A(n, m + n) = R(n)A(n, m) 
R{n) = S(n) - r(n)M(n) _1 iV(n). 
$(n) = 0(ra).H(n)0(7i)~ 1 , 

= 0(n)A(n, m + jn) 

= 0(n)P(n)v4(n, m + (j — l)n) 

= 0(n)P(n)0(n) _1 0(n),4(n, m + (j - l)n) 

= *(n)E/-(n) W( ._ 1)B 



(39) 
(40) 
(41) 



Wt/(n) 4m . (42) 



Remark: Note that the recursion operator (|4ip is nonlocal, but it does not generate the nonlocal 
higher flow equations, because in our cases, all the flow equations in (|26p are local and the recursion 
operator (|4ip is just extracted from these local flow equations. 

4. Applications 

In this section, we give some examples for the applications of the formula (|4ip for the recursion 
operator. Here we only consider 2-reduction and 3-reduction. 
2-REDUCTION 



For k = case, one calculates 



-ao(2) 

So the recursion operator [HKH] is 



a (2)=2n 2 , 0(2) = 8, 5(2) = a (2), T(2) = (0,l), 

M(2) = 1 - M j ' M(2) = ( \ 

( -d 2 \ 1 
N(2) =[ ] , i?(2) = -<9 2 + 2« 2 - cT W 



$(2) = -d^ua + ^d" 1 . (43) 
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U2t 3 = *(2)«2t a = 7«2xxx 4- 3u 2 u 2x , 



Beginning from «2ti = u 2x , we find 

1 

4 



15 2 5 5 1 

«2t 5 = 3>(^)«2t 3 = —U 2 U 2x + ~U 2 U 2xxx + -■U2XXW2X + — VQxxxxx- 

2 4 2 lb 

Note that in the reference [9], W.Strampp & W.Oevel also calculate the M(2) and N(2), but they are 
different from here. This is because we only use a,j{m). 
For k = 1 case, one has 



ai(2) = 2m, 0(2)=d, 5(2) = 0, r(2) = (ai(2), 1) 



29 \ , ( -%d 



a 1 (2) x -a 1 (2)-9 2 -29/ V i i( 2 ) 5_1 + ? / 

at(2) = ( - ai(2 o } " 92 ) , i?(2) = ia 2 - ^(2)3-^(2)0. 

The corresponding recursion operator [1] is 

<$>{2) = - A d 2 -u\-u lx d- 1 u 1 . (44) 



Thus, from u\ti = u\ x , we know 



uit 3 = $(2)ui tl = -u\ xxx - -u\ui x 



1 3 , 

15 a. 5 2 1 5 ^ 5 



M «5 = ®( 2 ) U U 3 = -5-«l«la: - o'Wirrrrrr'Wl + rrzU\ xxxxx - -U lx - -UiUi xx Ui 3 
o o Id o z 

At last, for k = 2, we have 
O(2) = n a-u 0l = ^ - 1 , 0(2)- x =u 9- 1 «o 2 > 5(2) = 0, T(2) = (0,n 2 ), 
M(2) = ' t 1 , M(2)" ] ' - 



-u, 



Id 2 -2u (u d + uqx) J V i n o 1(9 1-u o9m 2 ~lu 1 d 1 u 1 

10 



N(2) 



Wd 2 



1 



, R{2) = -u^u^d 3 . 
/ 4 



Therefore the recursion operator |31j is 



$(2) = -4d d u d^u 2 . 



(45) 



So 



«0ti 



0: 

*(2)«ot a 



1 



$(2)^0*3 = ■^UaO-O'U-OUOxxUOxxx + 5u 0xxx ul x + lOuoUOxxxx^Oa; + ^UqU 0xxxxx ). 



Obviously, all of above soliton equations are consistent with flow equations of (|19p ~ (|24p . which shows 
the validity of the explicit recursion operators (|43j) -([45|). 
3-REDUCTION 



k = case 

0(3) = 

r(3) = i 

M(3) _1 



d \ 
d J 

10 

oi(3) 1 

/ 



S(3) 



M(3) 



1 
3 



ao(3) - ai(3) x ai(3) \ 

ao(3)x + ai(3)xx a (3) - 2ai(3)a; / 

/ -3d 

-3d 2 -39 



\ 




\ -301(3)^-01(3)9 -d 3 -3d 2 -3d J 



3^ 
1 

3 



\ 



\ -|9 + ioi(3)9- 1 -|9- 1 ai(3) 

/ - ai (3) x - ai (3)d-d 3 -3d 2 \ 

N(3)= -00(3^ + 01(3)^ -201(3), -ai(3)9-9 3 

V 00(3)^-01(3)^ -200(3)^ + 301(3)^ / 

with ai(3) = 3u2, oq(3) = 3u^ + 3u2 X - Then by (f40l) and (liTj) . the recursion operator is given by 



d>(3) 



$11 $12 
$21 $22 



(46) 



where 



$11 
$12 
$21 



^da (3)d- 1 - IdatWxd- 1 + ~0ai(3) + ^ 3 + ~o (3) - ^1(3)*, 
^(S^ + ^^ + io^S), 

-^doo^cT 1 + ^oi^cT 1 - Hs3 0l(3)a -i _ ^4 _ 2 ai(3)aai(3)9 -i 

11 



2 11 

--ai(3)<9 2 - -00(3)3; + -ai(3) xx , 

$22 = ^da (3)8- 1 -Ida^d- 1 -^- |ai(3)5 + isai (3) + |a (3) - c^)*. 
We have checked the action of recursion operator (j46j) on the ti flow to £4 flow, that is, 



'«2 
'«3 



'«2 
U3 



U 2 
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^u 3 u 2x + ^u 2 u 3x + 2u 2 u 2xx + 2u 2 + \u 3xxx + ±u 2 



-2u 3x u 2x - 2u 2 u 3xx - 2u 2 u 2 



2x 1 3 

1 _ 2 

j'.rr 3^32:2:2:2: g U 2xxxxx 



Au\u 2x + 4u 3 u 3x - 4:U 2xx u 2c 



k = 1 case, 



0(3) = 
T(3) = I 

M(3) = 
M(3) _1 
iV(3) = 



9 ] 5(3)= f oi(3)-02(3) x \ 

v 5 y ' v -01(3),, +02(3)^ / ' 

a 2 (3) 1 \ 

, a 1 (3)-2a 2 (3) x a 2 (3) 1 J ' 

/ -3d 

-202(3)^- 202(3)9 -35 2 -3d 
\ -2ai(3) a + 302(3)^ - oi(3)fl - a 2 (3)<9 2 - d 3 -3a 2 (3) x - 2a 2 (3)d - 3d 2 

\ 



\ 


-39 j 



( -Id- 1 



A 

\ B 



3 a 

c 







/ - ai (3)d - a 2 (3)d 2 - d 3 



-a 2 (3) x -2a 2 {3)d-3d 2 



-a 1 (3) x + a 2 (3) xx - oi(3)S - o 2 (3)<9 2 - 9 3 
a i( 3 )xx -°2(3) 

XXX 



where 



ai(3) = 3u 2 + 3u\ + 3ui x , a 2 (3) = 3ui, 

A = Haa^a^ + i C = ia 2 (3)9- 1 -^- 1 a 2 (3) + i 

S = ^a 1 (3)a- 1 -i^V(3)-^a 2 (3)9- 1 -^a 2 (3)9- 1 + ^a-%(3)9a 2 (3)^ 
^a 2 (3)-^- 1 a 2 (3)a-^. 

12 



Then according to (|40p and (|4ip . we get the recursion operator 



/ $n $i 2 \ 

*(3) = " ' (47) 

V $21 $22 / 



where 



$n = -^a 2 (3)a- 1 a 1 (3)-^a 2 (3)a- 1 a 2 (3)a + ^a 2 (3)a + iaa 1 (3) + ^ 3 , 
y y y o o 

$12 = ^a 1 (3)a- 1 -^a 2 (3)^ 1 -^a 2 (3) 2 a- 1 -^a 2 (3)a- 1 a 2 (3) + ^ 2 , 

$21 = -~3ai(3)<9 - ^ai(3)0 2 - -flai^S"^!® - i a ?(3) - ^ai(3)a~ 1 a 2 (3)a 

-~ai(3)a 2 (3)d + ^ 2 a 2 (3)5- 1 a 1 (3) - -^ 2 a 2 (3)d + ^ 2 a 2 (3)^ 1 a 2 (3)9 
y 9 9 9 

+^a 2 (3)aa 2 (3)a- 1 a 1 (3) + ^a 2 (3)da 2 (3)d- l a 2 (3)d + A a2 (3)<9a 2 (3)d 

-io2(3)9ai(3) - -U(3)da 2 (3)d - -L(3)d 3 - ^ 2 ai (3) - V, 
y y y y y 

$22 = —dax^d- 1 + ^davffl^d- 1 - ^da 1 (3)a 2 (3)d- 1 - ^da 1 (3)d- 1 a 2 (3) 
-£0ai(3) - ^ai(3)a 2 (3) :l .a- 1 - ~a x (3)a 2 (3) - ioi(3)S + -^a^cT 1 
+^ 2 a 2 (3)<r 1 a 2 (3) - ~d 2 a 2 (3) + ^-a 2 (3)da 2 2 (3)d- 1 + ^a^da^d' 1 a 2 (3) 
-\a 2 {3) ai {3) x d- 1 + ^a 2 (3)a 2 (3) :E:E a- 1 - ~a|(3)d - ^a 2 (3)9 2 + ^^(3)^ 
+iai(3)5 - ^a 2 (3)5a 2 (3) 3: a- 1 - Va^)^" 1 - h\ 

With the recursion operator (|4"7|) . we can get the £4 flow from ii flow, 

12 4 o 2 8 

1*1*4 = ^ UlxXXX + 3 U2xxx + 2u ^U\ x - gUla^l + 2u 2 U\ xx + 4u 2 U 22 . + -MlMlxra + ~UlxUl xx , 

U2u = -^l^ixxUlx ~ 8uiu 2 u 2x - Auju 2 u lx - Auiu 2 ui xx - \u x u 2x u Xx - ^u 2 u\ x - ^u 2xxxx 

2 2 ^ 2 ^ 3 2 8 

--piUlxxxxx ~ <±U 2 U1 X ~ -U Xxxx U x - ~U x U 2x - 2u 2 U X xxx r.'Ul'U'lxxxx n Ui x U\ xxx 
9 o o o o 

2 10 2 2 4 o 

--^Uiu^xxx ~ 2u 2xx u lx - —u 2x ui xx - 2u 2 u 2xx - 2u lxx - 2u 2x - -Wfj.. 

k = 2 CciSG, 

0(3) = J • °( 3 )^ = P J _ a ._i -1 - 

V J \P 01 P 02 P 03 J 

13 



M(3) 



' -3a 3 (3)<9 - a 3x 
An 



3a 3 (3)<9 - 2a 3x 
D -1A 















3 u ° u 







3a 3 (3)<9 - 3a 3x J 
\ 



M(3) 



-i 



,4 




V 



/ 



iV(3) 



/ Dio D u ^ 

A>o 
V o j 



with 



a 2 (3) = 3uo(ni + u 0x ), a 3 (3) = Uq 



^12 = 02(3) - 03(3)2:, Pi3 = a 3 (3), P i = -02(3)2; + 03(3)2:2:, 

P12 = a 2 (3) - 2a 3 (3) x , P 13 = o 3 (3) , 

Doi = -02(3), + a 3 (3)^ - 2a 2 (3)d - 3a 3 (3)d 2 , 

D-ifi = 02(3)2-3; - a 3 (3) xxx - a2(3)9 2 - a 3 (3)<9 3 , 

= -2a 2 (3)2: + 3a 3 (3)^ - 2a 2 (3)<9 - 3a 3 (3)<9 2 , 
Dio = -a 2 (3)a 2 - a 3 (3)a 3 , D n = -2a 2 (3)9 - 3a 3 (3)9 2 , 
^oo = -a 2 (3)a 2 -a 3 (3)a 3 , 

A = —v^V-^Dnu^d- 1 ^*, C = ~Uq 3 a- 1 D_i,i«o 2 5 _1 «o 



Then the recursion operator can be got with the help of (1401) and (141 H . 




where 
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$21 = 



-u lx u x Pi 2 d 1 u q 2 Diquq8 1 u 2 + u lx P 13 ADi u d 1 u 2 

-^ui x u Q 1 Pi2d' 1 UQ 2 Diid' 1 ui x d' l u Q 2 + u lx Pi 3 AD u d~ l u lx d~ l UQ 2 

1 1 

~-^u 1x Uq 2 Pi 3 d~ 1 UQ D Q Qd" ui x d" Uq 2 + -u dP Q iUQ d^u^ 2 D w u d~ 1 u 2 

-u dPo2AD w uod~ 1 UQ 2 - UQdP^BDxQUQd^UQ 2 
+^u dPoiu 1 d- 1 u 2 D n d- 1 u lx d- 1 u 2 - uodPo2AD 11 d' 1 u lx d~ 1 u 2 

-UQdP Q3 BDi 1 d~ 1 ui x d~ 1 UQ 2 + ^u dP 02 u 2 d~ 1 u 1 D 00 d~ 1 ui x d~ 1 u Q 2 
-uodPosCDood^u^d^UQ 2 , 

$22 = -^uixUq 1 Pi2d~ 1 Uo 2 D 11 d~ 1 Uo 1 + u lx P i3 AD 11 d- 1 v^ 1 

1 1 

--ui x Uq 2 Pi 3 d~ Uq 1 D 00 d~ Uq 1 + -u <9-Poi^o 1 d~ 1 u 2 D n d~ 1 u 1 

-u Q dP 02 AD 11 d- l Uo 1 - uodPozBDud- 1 ^ 1 
+^u dP 2U 2 d~ 1 u 1 D 00 d~ 1 u 1 - uodPosCDood^UQ 1 . 

The t\ and £4 flows are as follows, 
Uoti = "lti = 0, 

u t 4 = -«o(-12'Ui a; Mf - 12ui x u 0x ui + 2uqUi xxx + 2u ul xx + u o u 0xxxx + 4-ui x uou 0:c:c + 6u xUqU1xx 

+2u u 0xxx u 1 + 4:U u 0xxx u 0x - 2ui x u 2 0x - u Qxx Uq x - 6u 0xx ul - GuoxxUoxUi), 
mu = -qUo(24ux x u 0x uou 0xx + 6u u 0xxx uiu 0x - 30ui x uiu u 0xx - 6u 0xx ul x - 72uq xx uq x u\ 

3 2 3 2 3 2 

+3u ui xxa ; I + ISmix^oUoxxx - 36uoxxMi + 33«oui xx « 03 . - 12ui x u 0x - Mui x u\u Qx 
+lQulu 0xxxx u 0x - 18u ui xx ul + 18u 0xx ulu 0xxx + 2Au^ui xxUq xx - 42uoxxU1Uq x 
-36uiu u 2 x + 18ui x UqUi xx - 6uiu ul xx - 12u uq xxx u\ + 22u u 0xxx u 2 ]x + §UqUq xxxx U\ 
-72u lx u\ + 12w 0x u ul xx - UAui x u 0x uj + 2ulu 0xxxxx + ^ului xxx u\ + 2Au 2 ui xxx uq x ) . 

Because of the complication of the corresponding calculations, we have only checked uot 4 by the 
recursion operator (|48p from the t\ flow. 

Remark: Prom the examples above, one easily find that though these recursion operators contain 
the nonlocal terms, the flows generated by them are local. And also the forms of the recursion operators 
are much more complicated when n and k tend to larger. 
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